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ABSTRACT 

We present far-infrared (57—196 /im) spectra of 21 protostars in the Orion molecular clouds. These 
■ were obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel 

Space observatory, as part of the Herschel Orion Protostar Survey (HOPS) program. We analyzed the 
emission lines from rotational transitions of CO, involving rotational quantum numbers in the range 
Jup = 14—46, using PACS spectra extracted within a projected distance of < 2000 AU centered on the 
protostar. The total luminosity of the CO lines observed with PACS {Leo) is found to increase with 
increasing protostellar luminosity (Ltoi)- However, no significant correlation is found between Lqo and 
evolutionary indicators or envelope properties of the protostars such as bolometric temperature, T^oi 
or envelope density. The CO rotational (excitation) temperature implied by the line ratios increases 
with increasing rotational quantum number J, and at least 3—4 rotational temperature components 
are required to fit the observed rotational diagram in the PACS wavelength range. The rotational 
temperature components are remarkably invariant between protostars and show no dependence on 
Lboi , Tboi or envelope density, implying that if the emitting gas is in local thermodynamic equillibrium, 
the CO emission must arise in multiple temperature components that remain independent of Lhoi over 
two orders of magnitudes. The observed CO emission can also be modeled as arising from a single 
temperature gas component or from a medium with a power-law temperature distribution; both of 
these require sub-thermally excited molecular gas at low densities (n(H2) ^ 10^ cm~^) and high 
temperatures (T > 2000 K). Our results suggest that the contribution from photodissociation regions 
(PDRs), produced along the envelope cavity walls from UV- heating, is unlikely to be the dominant 
component of the CO emission observed with PACS. Instead, the 'universality' of the rotational 
temperatures and the observed correlation between Leo and Li,oi can most easily be explained if the 
observed CO emission originates in shock- heated, hot (T > 2000 K), sub-thermally excited (n(H2) 
< 10^ cm~'^) molecular gas. Post-shock gas at these densities is more likely to be found within the 
outflow cavities along the molecular outflow or along the cavity walls at radii > several 100—1000 AU. 
Subject headings: circumstcllar matter — ISM: jets and outflows — molecular processes — stars: 
formation — techniques: spectroscopic 



1. INTRODUCTION 

The early evolution of protostars is driv en by the 
comp e tition between infall and outflow (iShu et al.l 
T987[ IStahler fc Pallal [20051 : IMcKee fc Ostrikeri 120071 : 
HartmannI l2009f ) . Energetic processes associated with 
mass accretion and ejection heat up the surrounding 
gas t o temperatures of several 100 K to several 1000 K 
(e.g.. Ivan Dishoeck et al.|[2009[ ). The infalling envelope 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA 



material first lands on the p rotostell ar disk before ac- 
creting onto the centra.1 star (iCassen fc Moos manI 119811 : 
iNeufeld fc HollenbachI 119941: iHartmanm i2009lV Uhra- 
violet radiation produced by accretion onto a proto- 
star can heat up th e circumstellar gas t o a few 100 K 
(|Spaans et al.l 119951 : Ivan Kempen et al.l l2009a[ ). Pro- 
tostars drive powerful bipolar jets and outflows into 
the surrounding medium, producing shocks, that com- 
press and heat up the envelope or ambient mate- 
rial to several 1000 K even at relatively l arge dis- 
tances (> 1000 AU) from the protostar (iBachilleii 
[TqM IBachiller fc TafaUal [19991 : IReinurth fc Ballvl 120011 : 
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IHollenbach fc McKed [TOSl IKaufman fc Neufeldl 119961) . 
While the lowest rotational transitions (J < 3) of CO 
have been widely used to trace cold (T < 30 K) molecu- 
lar gas associated with protostars, the far-infrared (FIR) 
CO lines (J >1 4) are better diagnostics o f warm and hot 
molecular gas (iWatson et al. 1980: Storey et alJ 119811 : 
iGiannini et all 120011: iNisini et all moW . These latter 
lines can in principle be used to characterize the density, 
temperature, molecula r abundance and spatial extent 
of the emitting region ("Watson et al'ligSS"; 'Draine et al.' 
1983; HoUcnbach & McKcc 1989; Kaufman & Neufeld 
19961 ). The spatial distribution and physical conditions 



in the emitting gas can provide vital clues about the heat- 
ing mechanisms and therefore also on the various ener- 
getic processes associated with the earliest phases of star 
formation. 

The first detection of the FIR CO lines from an astro- 
nomical source was toward O rion-KL, with the K uiyer 
Airborne Observatory (KAO: IWatson et al.lll980() . Be- 
cause of low sensitivity, the KAO could only observe 
the bright est regions around luminous and rnassive 
protostars ('St orev et all 119811 : iStacev et aLlll982L 119831 : 
iJaffc ct al. 198^ These early observations of very mas- 
sive protostars suggested that FIR CO emission arises in 
gas hea ted by non-dissociat i ve, magnetohydrod ynamic 
shocks (jWatson et all [l985l: iStorev et all Il981l) . FIR 
CO lines toward several low-mass protostars were later 
observed wit h the Infrared Space Observatory (ISO; 
e.g., Nisini ct al. 1997; Ceccarelli et al."1998;'Nisini et alj 
1999: Bcnedcttini et al.. ,2000.: Saraceno et al. .1999a.bl: 



Giaiinini efaLl 120011: " iNisini et all l2002t Ivan Dishoeckl 
20041 ). Modeling based on the ISO data found 
that the observed CO emission originates in molec- 
ular gas within a projected distance of < 1500 AU 
from the protostar, with temperatures in the range 
of 200—2000 K and II2 densities in the range of 
10^-10^ cm-3 (ICeccareUi et al.l[T998l: 



Benedettini et al .1120001: IGiannini et al.|[200 1: Nisin i et al.l 
20021) . These studies suggested slow {vs ^ 20 kms ^) 
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non-dissociative shocks as the main mechanism for the 
heating of the molecular gas. However, based on the sim- 
ilarity between the obse rved CO line rati os of protostars 
and Herbig Ae/Be stars iCeccarellil ()2000[ ) proposed that 
CO J„p=14— 21 lines from protostars probably originate 
in the dense Photodissociation Regions (PDRs) created 
by UV photons from the central objects. Lack of spatial 
information (beam size ^ 80") and low spectral reso- 
lution of the ISO-LWS, however, prevented any further 
progress in disentangling the various components con- 
tributing to the FIR CO emission. 

With the Photodetector Array Camera and Spectrom- 
eter (PACS; Poglitsch ct al. 20 ij) onboar d the Herschel 
Space Observatory (Pilbr att et al.l l2010[ ). it has now 
been possible to observe the rotational transitions from 
the entire CO ladder from J=13-12 to J=49-48 in 
the FIR wavelength range (j van Ke nipen et al .' '2010a"!9; 
Fich et al.l 120101 : iHerczeg et al. 20121: iGoicoe chca ct al. 
2012| ) at high spatial resolution (beam size ~ 10" at 63 



/xm) and unprecedented sensitivity. Modeling of the ab- 
solute and relative line fluxes of more than 30 CO ro- 
tational transitions accessible with PACS can be used 
to place strong constraints on the density and te mpera- 
ture of the emitting medium fe.g. lNeufeldll201^ . How- 
ever, the CO lines observed with PACS are spectrally un- 



resolved (e.g Ivan Kempen et all l2010al lbl: IHerczeg et al.l 
I2OI2D rendering it difficult to separate out the contribu- 
tion from various velocity components of the emitting 
gas to the total observed line flux. Spectrally resolved 
observations of the lowest— J CO lines {Jup S 8) to- 
wards low-mass protostars in the submillimeter wave- 
lengths have shown the presence of narrow (FWHM < 2 
km s-i) and relatively broad (FWHM ~ 10-30 km s^^) 
velocity components in the profiles of these lines. The 
narrow component was fou nd to dominate the total li ne 
flux for these transitions (Ivan Kempen et aTI l2009bllal : 
Ivan Dishoeck et al.l 120091: lYildiz et al.ll2012D . For lines 
with Jup < 4, these studies attributed the narrow com- 
ponent as due to emission from passively heated pro- 
tostellar envelope; for CO J = 6—5 and CO J = 7—6 
transitions the narrow component was found to be dom- 
inated by emission from UV-photon heated gas. Early 
modeling of the unresolved CO lines observed with Her- 
schel/PACS {Jup > 14) towards a few protostars sug- 
gested that the emission from UV-heated PDRs domi- 
nat ed the line flux even in th e PACS range up to Jup = 
20 (jvan Kempen et al]|2010b[ ) . These studies have shown 
that two separate components are required to explain the 
CO lines observed with PACS. Emission from PDRs, pro- 
duced along the envelope cavity walls due to UV-heating, 
was found to be the dominant component of the lower— J 
hues {Jup = 14-20) accessible to PACS. The higher- J 
CO lines {Jup ^ 25) were found to be dominated by the 
emissi on from the small scale shocks along the cavity 
walls (|van Kempen et aT1l2010bD . A more detailed mod- 
eling of the CO emission from three protostars wi th sim- 
ilar bolometric luminosities bv lVisser eTall (poll ) essen- 
tially confirmed this picture, except that in these models 
the PDR component can dominate emission even for the 
high— J CO lines, if the density of the envelope is suf- 
ficiently low. Thus, those authors tentatively suggested 
an evolutionary trend, in which the CO emission is dom- 
inated by shocks in the youngest source in their sample 
and by UV-heated gas in the oldest source with low enve- 
lope density (jVisser ct al. 2012). Herschel/PACS spectra 
of only a few sources have been studied so far, however. 
Analysis of a larger sample of protostars with a range 
in luminosity, evolutionary status, and envelope prop- 
erties is required before the above suggestions can be 
confirmed. In this paper, we carry out such a study. 

Here we present FIR (57—196 /im) spectra of 21 pro- 
tostars in the Orion molecular clouds observed with Her- 
schel/PACS. We obtained them as part of the Herschel 
Orion Protostar Surve y (HOPS), aii open t ime Key pro- 
grani with Herschel (iFischer et al.l 120101 : iStanke et al.l 
I20T0I : IMegeath et al.11201 2'). Our sample of protostars in 
Orion is the largest such sample at a common distance 
for which FIR spectra over the entire PACS wavelength 
range have been obtained and analyzed so far. Our anal- 
ysis begins with a search for possible correlations between 
the observed CO emission properties and protostellar lu- 
minosity and evolutionary status. We then model the 
observed CO fiuxes and line ratios to derive the excita- 
tion conditions of the warm and hot circumstellar gas 
and investigate the physical origin of CO emission from 
protostars. Our sample is described in § [2l Observa- 
tions and data reduction are described in § [3l In § IH 
we present the observed properties of the CO emission 
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Table 1 

Protostellar properties 



HOPS TD 




DEC ( 12000"! 




^ bol 














IX g 


10 


5 35 9.005 


-5 58 27.5 


3 


50 


6 


11 


5 35 13.409 


-5 57 58.1 


8 


59 


37 


30 


5 34 44.062 


-5 41 25.8 


4 


87 


6 


32 


5 34 35.450 


-5 39 59.1 


2 


60 


6 


56 


5 35 19.466 


-5 15 32.7 


18 


47 


185 


60 


5 35 23.328 


-5 12 3.06 


21 


60 


6 


68 


5 35 24.305 


-5 8 30.59 


5 


99 


19 


84 


5 35 26.570 


-5 3 55.12 


45 


93 


6 


85 


5 35 28.183 


-5 3 40.93 


13 


157 


4 


87 


5 35 23.472 


-5 1 28.70 


27 


39 


185 


91 


5 35 18.914 


-5 50.87 


5 


38 


370 


108 


5 35 27.074 


-5 10 0.37 


57 


41 


185 


182 


5 36 18.833 


-6 22 10.2 


62 


63 


19 


203 


5 36 22.838 


-6 46 6.20 


15 


45 


74 


288 


5 39 55.944 


-7 30 27.9 


103 


57 


37 


310 


5 42 27.677 


-1 20 1.00 


11 


55 


19 


329 


5 47 1.606 


17 58.88 


3 


77 


4 


343 


5 47 59.030 


35 32.86 


4 


88 


6 


368 


5 35 24.725 


-5 10 30.2 


55 


154 


4 


369 


5 35 26.971 


-5 10 17.1 


20 


42 


90 


370 


5 35 27.629 


-5 9 33.48 


217 


79 
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^ Fiducial envelope density, calculated at 1 AU in the limit of no rotation (see 
IKenvon et"an[T993) . 



from protostars in Orion and discuss various correlations 
found between CO emission and protostellar properties. 
In § [5j we model the observed CO rotational diagrams 
to constrain the physical conditions of the emitting gas 
and discuss the results. The physical origin of the FIR 
CO emission from protostars is explored in § |51 Finally, 
our conclusions are summarized in § [T] 

2. HOPS SPECTROSCOPIC SAMPLE 

The HOPS sample for the PACS spectroscopic ob- 
servations was selected from 300 Spitzer identified pro- 
tostars in Ori on with detectable 24 micron emission 
(jKrvukova et a l. 2012; Mcgcath ct al. 2012). We first se- 
lected objects with F^, > 7 mJy at 24 /im with Class 
or Class I spectral indices and other strong indications, 
such as deep silicate and ice features in the Spitzer IRS 
spectra, that they are indeed protostars. In this list 
we searched for objects that resemble, at infrared wave- 
lengths, the nearly face-on protostar NGC 1333 IRAS 4B: 
evidence of compact (< 10 " from the continuum posi- 
tion at > 24 /im) scattered light or H2 line emission at 
2.3, 3.6, 4.5 or 24 /im and no evidence for long jets or 
dust lanes suggestive of a near-edge-on orientation. Our 
motive here was for the sample to contain many objects 
for which our view to the dense inner regions would be 
relatively unextinguished. We initially identified 44 such 
protostars out of the parent sample of 300, all of which 
we observed with Spitzer IRS in its highest spectral res- 
olution (Watson et al. 2012, in preparation). From these 
44 sources, we chose 31 objects with bolometric luminos- 
ity uniformly sampled from 0.5 to 30 Lq to observe with 
PACS. An additional 5 objects with bolometric luminos- 
ity in the range of ^^20— 200 Lq were added to the PACS 
target list without regard to orientation. In total, 36 
protostars were observed with PACS in the spectroscopy 
mode as part of the HOPS program. Here we present 
PACS spectra of the brightest 21 sources in our sample. 



roughly uniformly sampling bolometric luminosity in the 
range of -1-200 Lq. The FIR spectra of the full HOPS 
sample will be presented in a later paper. Due to the 
difficulty in determining the precise orientation of the 
protostars, many of the sources which were suspected to 
be close to face-on may actually be observed at a more 
inclined orientation. A future paper will use modeling 
of the SEDs as well as HST 1.6 fim images to better 
constrain the inclinations for the HOPS protostars. 

The observed spectral energy distributions (SEDs) 
of the protostars in our sample are presented in Ap- 
pendix [X] Through several ancillary observing programs 
that complement HOPS, we obtained extensive wave- 
length coverage for the sources in our sample resulting in 
well-sampled SEDs. From these SEDs, we computed the 
bolometric luminosity (Lboi) and bolometric temperature 
(Tboi) of each source. Ltoi was calculated by integrating 
the SED over frequency. Since the SEDs are sampled 
at a finite number of frequencies, we interpolated over 
the intervening frequencies fol lowing the m i dpoint inter- 
polation method described in lEnoch et al.l ()2009|) . The 
flux upper limits were removed and we interpolated over 
them before integration. Further, the SEDs were ex- 
trapolated from the non-upper limit fiux at the longest 
observed wavelength (in most cases 870 /im or 350 /im) 
using (X u'^. Tboi, which is defined as the temperature 
of a blackbody with the same mean frequency as the 
source SED, was computed from the mean freque ncy of 
the source SED, following the m ethod of Mvcrs fc Laddl 
(l993) and lEnoch et al.l pOOOO . Given our wide wave- 
length coverage, and, in particular, given that the peak 
of the SED is well characterized, our integration method 
accounts for most of the luminosity, and the Lhoi and 
Tfioi are well constrained. The dominant source of uncer- 
tainity in Lhoi and Ti,oi is from the finite sampling errors, 
typically ~10— 15%. One of the sources in our sample, 
HOPS 369, has been modeled as a double source, i.e., a 
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Figure 1. Distribution of Lf,^; and Ti,^i of the p rotostars for whi c h PACS spectra a re presented. Ti,^i of 70 K is generally taken as the 
dividing line between Class and Class I sources l lChen et al.lll995l : lEvans et al.ll2009l) 
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Figure 2. pi (see text, § 2), as a function of Ltoi and Tfto;- 

disk-dominate d source p l us a c older, envelope-dominated 
source by Ad ams et alj ()2012l the source SOF 4 in the 
paper) . The envelope source dominate the observed 
SEP only at wav elengths > 37 /im (see Figure 2 in 
lAdams et al.ll2012l ): so we integrated the SED longward 
of 37 /xm to obtain Ltoi and Ttoi for HOPS 369. The Lboi 
and Tboi computed for the protostars in our sample are 
hsted in Table [H 

The distribution of the bolometric luminosities and 
temperatures of the protostars in our sample are pre- 
sented in Figure [T] Protostars in our sample span a 
large range in Lboi; from ~ 1.9 L© to 217 Lq. They 
generally are in different evolutionary phases. Protostars 
are generally classified in to Class and Class I sources, 
two different, but related, evolutionary phases. Class 
sources are the youngest protostars and have most of 
their mass still in the infalling envelope (|Andre et al.l 
119931 I Andre fc Montmerld [1991 : they are highly ob- 
scured objects, invisible - in some cases - even at mid- 
infrared wavelengths. The more-evolved Class I sources 
have envelopes of mass less than the protostar itself 
and are visible in t he near-infrared wavelengths (jLadal 
[iMTt iWilkintd [19891 ). Tboi of 70 K is generally taken as 
the dividing line between Class and Class I sources 
(|Chen et al.lll995l: lEvans et al.ll2009[ ). Based on this cri- 
terion, our sample has 13 Class sources and 8 Class I 



sources. 

To follow the evolution of the physical conditions in 
protostellar envelopes, we estimated the envelope densi- 
ties of the protostars in our sample by modeling their 
SEDs. We assu me the density f ollows the rotating col- 
lapse solution of iTerebev et al.l ()1984[ ). While the den- 
sity varie s throughout t he en velope, we adopt the ap- 
proach of lKenvon et al.l ()1993l ) and use a fiducial density 
pi, calculated at 1 AU in the limit of no rotation, to 
track the envelope density. We computed pi by fitting 
the observed SEDs of the protostars with a large grid of 
models created with the Mon te Carlo radiative transfer 
code of I Whitney et al.l (|2003[ ). The detailed description 
of the grid and th e fitting procedu re ca n be found in 
iTobin et alj (|2008| ). lAli et al.l (|2010f ) and iFischer et all 
f|2010., 20121 ). For six sources in our sample (HOPS 60, 
68, 108, 368, 369, 370) SED model i ng re s ults have al- 
ready been published (jAdams et al.l 120121 : iPoteet et al.l 
1201 If ): for these objects, values of pi obtained from 
our modeling agreed with those published within a fac- 
tor of 2-3. For the composite SED of HOPS 369, we 
adopted the Oi valu e derived for the envelope source by 
lAdams et al.l ()2012|) . The envelope pi values for our sam- 
ple are also listed in Table [1] and are plotted as a func- 
tion of Lboi and Ttoi in Figure [2j As can be readily seen, 
pi appears to be uncorrelated with Ltoi- The Spear- 
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Table 2 

Log of PACS spectroscopy observations 



HOPS ID 


OBSID 


OD 


Date 


Total time 
(s) 
(s) 


Observing mode 


Primary wavelength 
ranges observed 
(/im) 


10 


1342215693 


665 


10 Mar 2011 


3079 


Pointed /Unchopped 


57-71 


& 


102-142 




1342215694 


665 


10 Mar 2011 


1930 


Pointed /Unchopped 


71-98 






11 


1342204115 


482 


07 Sep 2010 


5156 


Pointcd/Chop-Nod 


71-98 


& 


102-142 




1342204116 


482 


07 Sep 2010 


3369 


Pointcd/Chop-Nod 


57-71 






30 


1342215697 


665 


10 Mar 2011 


3079 


Pointed/Unchopped 


57-71 


& 


102-142 




1342215698 


665 


10 Mar 2011 


1930 


Pointed /Unchopped 


71-98 






32 


1342192115 


303 


12 Mar 2010 


3129 


Poiiitcd/Chop-Nod 


57-71 








1342192116 


303 


12 Mar 2010 


5020 


Poiiitcd/Chop-Nod 


71-98 


& 


102-142 


56 


1342227334 


833 


24 Aug 2011 


3079 


Pointed/Unchopped 


57-71 


& 


102-142 




1342227335 


833 


24 Aug 2011 


1930 


Pointed /Unchopped 


71-98 






60 


1342227611 


839 


30 Aug 2011 


8933 


Mapping/Unchopped 


57-71 








1342227613 


839 


30 Aug 2011 


8933 


Mapping/Unchopped 


57-71 








1342227616 


839 


30 Aug 2011 


6083 


Mapping/Unchopped 


71-98 


& 


102-142 




1342227618 


839 


30 Aug 2011 


6083 


Mapping/Unchopped 


71-98 


& 


102-142 


68 


1342226202 


823 


14 Aug 2011 


3079 


Pointed/ IJjjcJioppcd 


57-71 


& 


102-142 




1342226203 


823 


14 Aug 2011 


1930 


Pointed/Unchopped 


71-98 






84 


1342227338 


833 


25 Aug 2011 


3079 


Pointed/Unchopped 


57-71 


& 


102-142 




1342227339 


833 


25 Aug 2011 


1930 


Pointed /Unchopped 


71-98 






85 


1342215657 


664 


8 Mar 2011 


8933 


Mapping/Unchopped 


57-71 








1342215659 


664 


8 Mar 2011 


8933 


Mapping/Unciioppcd 


57-71 








1342215662 


664 


8 Mar 2011 


6083 


Mapping/Unchopped 


71-98 


& 


102-142 




1342215664 


664 


9 Mar 2011 


6083 


Mapping/Unchopped 


71-98 


& 


102-142 


87 


1342215650 


664 


8 Mar 2011 


3079 


Poiiitcd/Uiichoppcd 


57-71 


& 


102-142 




1342215651 


664 


8 Mar 2011 


1930 


Pointed/Unchopped 


71-98 






91 


1342215654 


664 


8 Mar 2011 


3079 


Pointed /Unchopped 


57-71 


& 


102-142 




1342215655 


664 


8 Mar 2011 


1930 


Pointed /Unchopped 


71-98 






108 


1342239690 


1020 


27 Feb 2012 


3079 


Pointed /Unchopped 


57-71 


& 


102-142 




1342239693 


1020 


27 Feb 2012 


1930 


Pointed /Unchopped 


71-98 






182 


1342226753 


826 


18 Aug 2011 


3079 


Pointed/Unchopped 


57-71 


& 


102-142 




1342226754 


826 


18 Aug 2011 


1930 


Pointed/Unchopped 


71-98 






203 


1342191363 


290 


28 Feb 2010 


3129 


Pointed/Chop-Nod 


57-71 








1342191364 


290 


28 Feb 2010 


5020 


Pointed/Chop-Nod 


71-98 


& 


102-142 


288 


1342227754 


838 


29 Aug 2011 


8933 


Mapping/Unchopped 


57-71 








1342227756 


838 


29 Aug 2011 


8933 


Mapping/Unchopped 


57-71 








1342227759 


838 


29 Aug 2011 


6083 
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Figure 3. PACS spectral map of HOPS 182 in tlie B2B 
(72—98 fim) spectral band. Within the square boxes, which rep- 
resents the 5x5 spaxel array, the observed spectra at that loca- 
tion is shown. The PACS spectroscopic field is overlaid on the 
PACS 70^01 image of HOPS 182 obtained as part of the HOPS 
program. The star symbol represents the MIPS 24/im coordinates 
of HOPS 182, which was used for telescope pointing. 

man rank coefficient, r^, for the correlation is 0.27 and 
the probabihty that that this correlation could be pro- 
duced by two random variables of the same sample size, 
P{rs), is 23%. On the other hand, pi is tightly cor- 
related with Tboi] the Spearman correlation coefficient, 
rs = - 0.82 and P{r,) < 0.001%, indicating that the 
correlation is statistically significant. The envelope den- 
sity generally decreases with increasing Ti,oi suggesting 
that Tboi is a good indicator of the evo lutionary stage of 
the protostars in our sample (also see lChen et al.|[l995l : 
lEvans et al.|[2009l) . 

3. OBSERVATIONS AND DATA REDUCTION 

The log of Herschel PACS observations is provided 
in Table [21 All spectra were obtained in the range 
scan mode to achieve a total wavelength coverage of 
57 /im to 196 /im. The wavelength range of 57—71 /.tm 
was observed in the 3'"'' grating order (B3A), the range 
71-98 ^m in the 2"^^ order (B2B), and the range 
102—196 iim in the 1** order (Rl). Nyquist sampling, 
with grating step size of 6.25 spectral pixels was used to 
cover the broad wavelength ranges in each spectral band. 
Seventeen of the sources were observed in the pointed 
mode and four were observed in the full PACS spatial 
resolution mapping mode. The pointed observations of 
the bright sources {F^q 10 Jy) were carried out in 
the unchopped mode and those of the fainter sources 
{Pio ^ 10 Jy) in the chop/nod mode. The mapping 
observations were done in the unchopped mode. For the 
chop/nod observations, a maximum chopper throw of 6' 
was used. For unchopped observations, off-positions free 
of cloud emission within 2° from the target were used 
to subtract the instrument background. These emission- 
free off-positions were determined from the CO J = 3 — 2 
maps of the region (Di Francesco et al. 2012, in prepa- 
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Figure 4. Continuum subtracted PACS spectra of 8 protostars 
which are representative of the entire sample. The central wave- 
lengths of the observed CO lines ( Jiip=14— 46) are indicated at the 
top. The spectra have been offset for better viewing. The contin- 
uum has been subtracted by applying a median filter to the spectra. 
An arbitrary offset has been added to the spectra to separate them 
in the plot. 

ration). 

All spectra were reduced using the Herschel Interactiv e 
Processing Environment (HIPE) version 8.0 f otMmxh . 
We used the pipeline scripts to obtain the final spectral 
maps over an area of 47" x 47" centered on each source. 
An example of the PACS spectral map is shown in Fig- 
ure [21 where the 5x5 spatial pixel ( "spaxel" ) array and 
the spectrum corresponding to each spaxel is displayed. 
The telescope pointing in all our observations was good: 
we mapped the continuum in each of the PACS spec- 
tral bands and found that the continuum peak was well 
centered on the central spaxel within the pointing uncer- 
tainty of ~ 2". The 1-D spectrum used in the analysis 
was extracted from the central spaxel (9.4" x 9.4") where 
the continuum emission from the protostar peaked. For 
the mapping observations, we extracted the 1-D spec- 
trum within an 9.4" x 9.4" aperture centered on the con- 
tinuum peak. The FWHM of the spectrometer beam is 
a function of wavelength and ranges from 10" at 75 /xm 
to ~ 12" at 160 ^J.m (see PACS Observer's Manual). We 
applied a wavelength dependent PSF-loss correction to 
the extracted spectrum to account for the different frac- 
tions of the PSF structure seen by the central spaxel at 
different wavelengths. 

Representative FIR spectra (57—196 /im) of a few 
sources in our sample are shown in Figure [4l At the 
distance of Orion (^ 420 pc), the extracted 1-D spectra 
corresponds to the emission from within a radius of ^ 
2000 AU from the protostars. Our PACS spectroscopic 
observations were designed to cover the rotational transi- 
tions in the ground vibrational state of CO ranging from 
J = 14 — 13 to J = 46 — 45. Several of these lines are 
seen in emission in the final extracted spectrum of most 
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Figure 5. Measured CO line widths (left) and l-cr line sensitivity of our observations (right). Lines detected at 3-cr level or higher are 
shown. Pointed/unchopped (black solid circles), pointed/chop-nod (blue circles), and mapping/unchopped (red circles) observations are 
identified. The gray solid lines correspond to HSpot predictions for our observational set up. The green solid lines are fits to the observed 
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Figure 6. Optical depth, t, of the CO J = 14 — 13 line as a 
function of observed line fiux for three angular sizes (for d = 420 
pc to Orion) of the emitting region. The optical depths shown are 
for a CO Une width {FWHM) of 10 kms"!. 

objects in our sample, indicating the presence of warm 
CO gas associated with these protostars. In four ob- 
jects, HOPS 11, 30, 91, and 329, we do not detect any 
CO lines above the 3a rms level in the PACS wavelength 
range. In the remaining 17 objects, several CO lines are 
detected, ranging from 3 lines in HOPS 84 up to 30 lines 
in HOPS 370 and HOPS 108 (see Figure gl). We mea- 
sured the fluxes of the observed CO lines in the spectra 
of the protostars by fitting them with Gaussian profiles. 
The measured CO line widths and the line sensitivities 
achieved from our observations are shown in Figure O 

4. OBSERVED PROPERTIES OF FIR CO 
EMISSION FROM PROTOSTARS 

4.1. CO excitation diagrams & rotational temperatures 

Since several CO lines are detected for many objects 
in the PACS wavelength range, it is useful to construct 
a rotational excitation diagram, which shows the rela- 
tive populations of CO rotational states derived from 
the observed line fluxes. In this diagram, the natural 
logarithm of the detected number of CO molecules in a 
rotational state divided by the degeneracy of that state, 
\n{Nj/gj), is plotted as a function of the energy, Ej (see, 
e.g., [Goldsmith fc Langer 1999, ). If the CO emission is 



optically thin then the number of CO molecules in the 
J*'' rotational state is given by, 

Nj = {4 7r cP Fjj^,)/{h Kj.j-i 

where Fj^j-i, vjj-i and Aj^^i are the line flux, fre- 
quency and Einstein A-coefficient, respectively, corre- 
sponding to the transition J — > J—1 and d is the distance 
to the source. The FIR CO lines observed with PACS 
are likely to be optically thin if the size of the emitting 
region is sufficiently large (^1")- Figure[6]shows the opti- 
cal depths of the CO J = 14—13 line observed toward the 
protostars in our sample for three assumed angular sizes 
of emission. To compute the optical depth, we assumed 
an excitation temperature of 300 K and a line width 
[FWHM) of 10 kms-i. The observed CO J = 14 - 13 
transition, which is the lowest excitation line observed 
and the most likely one to be optically thick, has r 1 
only if the size of the emitting region is <C 1" (<C 420 AU 
at the distance of Orion). The extent of the observed CO 
emis sion is likely to be much larger (~ 1000 AU) than 
this ("Viss er et all 120 121) . Moreover. lYildiz et all ()2012f ) 
have shown that even for the CO J = 6—5 line the op- 
tical depths are low (r < 2) for broad velocity compo- 
nent {FWHM > 10-15 km s~i). The higher excitation 
{Jup > 14) CO lines observed with PACS are likely to 
be broader than this, which would make the line optical 
depths even lower. Thus, the optically thin assumption 
appears valid and the above equation can be used to com- 
pute Nj and to construct the rotational diagrams. The 
CO rotational diagrams for the 17 protostars with CO 
detections are shown in Figure [Tj 

The observed rotational diagrams of the protostars 
typically show a positive curvature, which is more 
apparent in cases where more than 10 lines are de- 
tected. The rotational diagrams of a few other 
protostars whose Herschel/VKCS sp ectra have been 
published (Ivan Kempen et al.l l2010bl: iFich et all 120101: 
iHerczeg et al.l 120121 ) also sho w simi lar positive cur- 
vature as noted by iNeufeldl (|2012f ). This means 
that the rotational temperature defined as Trot — 
— {k d\n[Nj/gj\/dEj)~'^ increases monotonically with 
Ej. An average rotational temperature can be com- 
puted for a relatively small range in Ej by fitting 
a straight line to the points in the observed rota- 
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Figure 7. CO rotational excitation diagrams for the 17 protostars towards which CO lines are detected. Linear fits to the CO fines in the 
Rl (550 K < Ej/k < 1800 K; Jup = 14-25) band (red solid fine), the B2B (2000 K < Ej/k < 3700 K; Jup = 27-36) band (green sofid 
line), and the B3A (3700 K < Ej/k < 6000 K; Jup = 37—46) band (blue solid line) and the average rotational temperatures derived, 
Tri, Tb2B, and Tb3A are shown. Downward arrows indicate upper limits to the fluxes of CO lines which are blended with a nearby line. 
Open circles correspond to 3cr upper limits for the non-detections. Sources in which no CO lines are detected are not shown. 



tional diagram. We begin our analysis by comput- 
ing the average slopes of the observed rotational dia- 
grams for the PACS spectral bands Rl (102—190 /im; 
550 K < Ej/k < 1800 K; J„p = 14-25), B2B (72-98 
/im; 2000 K < E,/k < 3700 K; J^p = 27-36), and B3A 
(57-71 Aim; 3700 K < Ej/k < 6000 K; J„p = 37-46). 
The best linear fit for the average slope was obtained by 
minimizing x^l from these slopes, we obtained the aver- 
age rotational temperatures, Tri, Tb2b and Tbsa which 
are shown in Figure [7] and are listed in Table [3l This 
method of computing Trot for CO lines observed within 



the same spectral band has the advantage of minimizing 
possible systematic uncertainties introduced by any mis- 
match in the relative fiux calibration between different 
spectral bands. 

Figure [8] shows the distribution of the rotational tem- 
peratures Tri, Tb2b, and Tbsa- 7ri has a narrow dis- 
tribution ranging from 144 K to 375 K with a median 
value of 294 K. The uncertainties in Tri are relatively 
small with an average value of 21 K. Tb2b ranges from 
550 K to 924 K and has a median value of 739 K; it is 
less well determined and has relatively large uncertain- 
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Figure 7. continued 
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Table 3 

CO Rotational temperatures and luminosity 





Tlri 




Tri 


7b2B 


Tb3A 


Leo 


HOPS ID 


(K) 


(K) 


(K) 


(K) 


(K) 


(X 10-3 Lq) 



10 


238 ± 44 


391 ± 280 


315 ± 43 






1.3 


11 












< 0.7 


30 












< 0.8 


32 


307 ± 46 


282 ± 91 


275 ± 19 






1.8 


56 


218 ± 10 


338 ± 58 


262 ± 11 


550 ± 361 




8.2 


60 


221 ± 11 


303 ± 27 


294 ± 13 






6.2 


68 


144 ± 11 




144 ± 11 






1.6 


84 


186 ± 49 




186 ± 49 






0.7 


85 


186 ± 37 




270 ± 34 






1.2 


87 


270 ± 19 


307 ± 49 


316 ± 13 


575 ± 475 




7.6 


91 












< 0.8 


108 


290 ± 3 


410 ± 4 


344 ± 12 


734 ± 30 


1076 ± 25 


238.0 


182 


230 ± 3 


413 ± 15 


331 ± 20 


703 ± 27 


1220 ± 223 


48.4 


203 


322 ± 34 


443 ± 35 


350 ± 16 


773 ± 96 




9.7 


288 


234 ± 11 


362 ± 45 


315 lb 19 






10.0 


310 


252 ± 11 


375 ± 49 


268 ± 6 


924 ± 152 




9.8 


329 












< 0.7 


343 


260 ± 51 


570 ± 417 


375 ± 43 






2.0 


368 


193 ± 8 


386 ± 86 


262 ± 17 


739 ± 313 




6.4 


369 


191 ± 6 


398 ± 52 


264 ± 18 


807 ± 168 




11.0 


370 


239 ± 4 


374 ± 17 


332 ± 13 


807 ± 42 


1775 ± 271 


58.0 



Note. — Tlri is computed for the CO lines in the range J„p=14— 18, Tg^i for 
J„p=19-25, Tri for J„p=14-25, Tb2B for J„p=27-36 and Trsa for J„p=37-46. 
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Figure 9. Distribution of rotational temperatures Tlri (red solid 
bars) and Tsm (blue line bars). The dashed lines represent the 
median values of Tlri (red) and TgRi (blue). 

ties (185 K). Tbsa could be determined only for three 
sources and ranges from 1076 K to 1775 K (see Table [3] 
and Figure[5]). Since, as pointed out above, Trot increases 
monotonically with Ej (or J), the estimated average ro- 
tational temperatures generally increase with the number 
of lines detected in a spectral band. The spread in rota- 
tional temperatures listed in Table [3] is partly due to this 
behavior. For instance, for eleven sources for which more 
than 9 CO lines are detected in the Rl band, the corre- 
sponding Tri values show a much smaller range (262 K 
to 350 K with a median value of 315 K). Similarly, for six 
sources in which 5 or more lines are detected in the B2B 
band, the corresponding Tb2b values ranges from 703 K 
to 924 K. Thus, the observed CO rotational tempera- 
tures, in particular Tri, is very similar for the protostars 
in our sample. 

We also computed the average rotational tempera- 
tures separately for the PACS spectral bands Long Rl 
(142-190 fim- 550 K < Ej/k < 1000 K; J„p 14-18) 
and Short Rl (103-140 ^m; 1000 K < Ej/k < 2000 K; 
J„p = 19—25) which we will call Tlri and TgRi. The 
rotational temperatures are computed only if there are 
more than two lines detected in a given spectral band. 
The distribution of Tlri and TgRi for sources where both 
these temperatures could be determined are presented in 



Figure IHl Although there is some overlap between the 
derived values of Tlri and Tsri, for most objects there 
is a significant difference between the two. The median 
value of Tlri is 239 K and that for TgRi is 386 K; this 
difference is more than that can be accounted for by the 
uncertainties in these quantities. This difference suggests 
that even in the Rl spectral band there is a significant 
curvature in the rotational diagram for most objects in- 
dicating that even the mid— J CO lines (J„p = 14 to Jup 
= 25) observed in the PACS Rl spectral band cannot 
be explained by a single temperature component in local 
thermodynamic equilibrium (LTE). 

4.2. Correlation with protostellar properties 

The CO emission lines detected in the PACS wave- 
length range have their origin in the warm gas in the im- 
mediate vicinity (i.e., at projected distance < 2000 AU 
at 420 pc) of the protostars. Since the energy that goes 
into the heating of the circumstellar gas has to be derived 
ultimately from the central protostar, it is instructive 
to search for possible correlations between the observed 
characteristics of the CO emission and protostellar prop- 
erties. 

We first computed the total observed CO luminosity, 
-^CO, by adding up the fluxes of all the lines detected 
with S/N > 3a. Lqo is plotted as a function of Ltoi, 
Tboi, and the envelope density at 1 AU, pi, for all the 
protostars in our sample in Figure 1101 For sources where 
no CO lines were detected within the PACS wavelength 
range, an upper limit to Tco was computed by adding 
up the 3ct detection limit at the location of each CO line. 
To carry out a correlation analysis accomodating the up- 
per limits, we use d the Astronomical Survival Analysis 
(ASURV) package (Lavallev et al."1992) that implements 
the method s prese nted in Feigelson & Nelson (1985) and 
llsobe et al.l (|1986l ). We used the generalized Kendall's 
tau and Spearman's rank order tests to compute the cor- 
relation probabilities between Tco and the protostellar 
properties. These tests show that ico is tightly cor- 
related with Lboi- The associated Spearman rank cor- 
relation coefficient, = 0.71 and the probability that 
Leo and Ltoi are uncorrelated is 0.2%; the generalized 
Kendall's tau test gives this probability to be 0.3%. On 
the other hand, no statistically significant correlation is 
found between Tco and the protostellar evolutionary in- 
dicators Tboi (Spearman prob. = 17% and Kendall's tau 
prob. — 18%) or pi (Spearman prob. = 8% and Kendall's 
tau prob. = 9%). 

CO rotational temperatures Tri and Tb2b (see Sec- 
tion UT]) are plotted as a function of Lboi, T^oi and pi for 
the protostars in our sample in Figure 1111 Tri is uncor- 
related with Lboi and remains constant over more than 
two orders of magnitude in Tboi- Tb2b could be deter- 
mined only for protostars with Lboi ^ lOT© as high— J 
CO lines (Ej > 2000 K) are detected only for these ob- 
jects. Tb2b also appears to be uncorrelated with Liioi and 
within the uncertainties it remains roughly constant over 
a smaller range in TboZ- Further, Tri and Tb2b are found 
to be uncorrelated with the Tboj and pi values of the pro- 
tostars, as shown in Figure [TTJ Additionally, Figure [12] 
shows that Tlri and Tsri also are uncorrelated with pro- 
tostellar luminosity; they remain more or less constant 
over a large range in Lboi- This behavior demonstrates 
that irrespective of the range in Ej (or J) over which 
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the average rotational temperatures are computed, they 
remain independent of Lboi- 

In summary, the observed CO hne fluxes and luminosi- 
ties scale with Ltoi- The relative line fluxes (as quantified 
by the rotational temperatures), however, appear inde- 
pendent of Lboi- Both the CO line luminosities and the 
CO rotational temperatures show no significant correla- 
tion with evolutionary indicators or envelope properties 
such as Tboi, and pi. 

4.3. CO emission from an outflow lobe 

In the PACS field observed toward one of the sources 
in our sample, HOPS 87, we detected CO emission from 
the extended lobe of a molecular outfiow away from the 
exciting source. The PACS spectral map centered on 
HOPS 87 is shown in Figure fT3l Anoth er protostar, 
HOPS 88 (a.k.a. MMS 5: IChini et all[l997l) . seen toward 
the north-western edge of this field, is driving a bipo- 
lar outflow along the east-w est direction (William s et al.l 
I2003t iTakahashi et al.ll2008D . The red lobe of this outfiow 
falls on the two spaxels east of the on-source spaxel where 
HOPS 88 is located. Figure [T3l shows that the line emis- 
sion is brighter in the outfiow spaxels labeled "off-source 
1 & 2" than in the on-source spaxel. The CO rotational 
diagrams for the on-source spaxel and both off-source 
spaxels corresponding to the redshifted outfiow lobe are 
shown in Figure The rotational temperature, Tri , 
measured for the on-source spaxel and for the outflow 
spaxels are nearly the same, within the uncertainties. 
The luminosities of the CO lines detected with PACS at 
the outflow lobe positions far away from the protostar 
are higher than that from the spaxel where HOPS 88 
is located. For example, Leo from the off-source 2 po- 
sition, which is ~ 19" (projected distance ~ 8000 AU) 
from the protostar, is about a factor of 2 higher than that 
for the on-source position (see Figure [H)) . Additionally, 
the rotational temperature, Tri, and the CO luminosity, 
Leo J observed for the emission from the off-source po- 
sitions are very similar to those found for the on-source 
positions of all the other protostars in our sample (see 
Figure El &[IOl). 

5. EXCITATION CONDITION OF THE CO 
EMITTING GAS 

Our most striking result is the independence of CO line 
ratios (or equivalently rotational temperatures) with re- 



spect to protostellar luminosity, and in a relatively large 
sample of protostars. The interpretation of this result 
and other apparent correlations, or lack thereof, between 
the observed properties of the CO emission and the pro- 
tostellar properties depends on the physical conditions 
(temperature and density) in the emitting molecular gas. 
The CO emitting gas will be in LTE if the total gas 
(molecular hydrogen) density, 7t.(H2), is greater than ricr, 
the critical density of the CO transitions at a given tem- 
perature. The critical density of the lowest— J transition 
( J = 14 - 13) that we detect is 2-3 x 10^ cm"^ and 
that of the highest— J transition ( J = 46 — 45) is 5—7 x 
10'' cm~'^ for t emperatures in t he range of 300—3000 K 
(|Neufeldll20ll lYang et aLllMol ). Therefore, at densities 
n(H2) > 10® cm~'^, the CO rotational states will be ther- 
malized and the observed rotational temperature is the 
physical temperature of the molecular gas. For gas densi- 
ties < 10^ cm~^, the CO excitation is sub-thermal\^ and 
the rotational temperature can be significantly different 
from the physical temperature of the gas. 

5.1. Modeling of CO emission from protostars 

Recently, iNeufeldl ()2012[ ) has demonstrated that the 
observed curvature in the CO rotational diagrams ob- 
tained with Herschel/PACS can be described surpris- 
ingly accurately by optically thin FIR CO emission orig- 
inating in an isothermal medium or a medium with a 
power-law distribution of temperatures, both at uniform 
density. In the following we use these sets of models to 
explore the likely physical conditions of the CO emitting 
gas in individual sources. 

5.1.1. Isothermal medium 

We first consider the simple case of an optically thin 
medium of uniform temperature and density. We com- 
puted the synthetic CO rotational diagrams for such a 
medium for a large range of values of T, n(H2), and a col- 
umn density parameter iV(CO) (for details see INeufeldl 
12012 ). We then compared the model rotational diagrams 
with the observed ones to find the best fit solutions by 
minimizing the reduced-x^. 

The best fit models and physical parameters (corre- 
sponding to the minimum in reduced-x^) obtained for 
CO emitting gas for all the protostars are shown in Fig- 
ure 1151 We only modeled sixteen sources for which more 
than three CO lines are detected. The best fit values 
of T and n(H2) are shown as a function of Lhoi in Fig- 
ure [TH For 5 sources for which less than 8 CO lines were 
detected, the density and temperature of the emitting 
medium are not constrained by the models. For all the 
other protostars, these parameters are well constrained; 
they are listed in TableH) The optimal solutions have gas 
temperatures T > 2000 K and densities ^(Hs) < 10^-^ 
cm^"^. While gas temperatures are well constrained for 

^ In tills paper, as in previous papers In the literature, tiie ex- 
citation is said to be sub-thermal If tiie level populations of the 
upper states are smaller than the values that would be obtained In 
LTE; conversely, the excitation Is described as thermal if the rela- 
tive level populations are given by Boltzmann factors appropriate 
to the gas temperature. In either case, the process of coUisional 
excitation Involves colliding molecules that are assumed to have 
a thermal (I.e. Maxwrell-Boltzmann) distribution of translational 
energies. 
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Figure 13. PACS field of the pointed observations centered on HOPS 87 overlaid on Spitzer/TRAC 4.5 fira map. The protostar HOPS 88 
(MMS 5: IChi ni et al.lll9971) is seen to the north-west within the spaxel labeled on-source. HOPS 88 is driving an outflow in the east-west 
direction. Bright H2 emission from outflow lobes is seen on either side of HOPS 88 in the IRAC image. The redshifted lobe of the outflow 
falls in the two spaxcls to the east, which are labeled as off-sourcel and off-source2. A continuum subtracted spectrum in the spectral 
band Long Rl (142—190 fim) is displayed in each spaxel. 
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Figur e 14 . CO rotational diagrams for the on-source spaxel and the off-source outflow spaxels associated with HOPS 88 as shown in 
Figure [T3l 
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Figure 15. (Left) Rotational diagrams of the observed CO emission with the best fits for the isothermal medium model overlaid (solid red 
line). Downward arrows indicate upper limits to the fluxes of the CO lines which are blended with a nearby line. Open circles correspond 
to 3(7 upper limits for the non-detections. (Right) Reduced-x^ contours for T and n(H2). The star symbol marks the minimum value of 
the reduced-x^. The contours corresponding to 68.3% (blue), 95.4% (green) and 99.7% (red) confidence levels are shown. 
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Table 4 

Best fit model parameters 



Isothermal Power law 



HOPS ID T n(H2) ""Leo (total) power law (i(H2) ""Leo (total) 

(K) (x 1(>' cm ■•) (x 10 L,.,) index, /) (x K)-"' em ■') (X 10 Lq) 
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Leo listed here is the total CO luminosity derived from best fit models for rotational transitions 
ranging from J = 1 — up to J = 80 — 79. 
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Figure 16. Best fit values of T and n(H2) for CO emission 
from an isotliermal medium plotted against Ltol ■ The uncertainties 
shown correspond to 99.7% confidence interval. The gray shaded 
regions indicate T & n(H2) parameter space not explored by our 
models. 



these sources, for many of them only an upper hmit could 
be obtained for densities. Also, several of them have the 
best fit densities well below 10'' cm"'^. This is because 
at high temperatures (T > 2000-3000 K) and low den- 
sities (n(H2) < lO'*'^ cm"'^), the CO line ratios (or rota- 
tional diagrams) are insensitive to density; at these high 
gas temperatures the rotational diagra ms are indistin - 
guishable for densities below 10 ° cm^ -3 (Neufeld'20ll)- 
Thus, our modeling shows that if the observed CO emis- 
sion originates from a single isothermal component, then 
the CO excitation is sub-thermal. 

For sources where the density and temperature of the 
CO emitting medium are well constrained, the total CO 
luminosities (for transitions J = 1 — up to J = 80 — 79) 
obtained from the best fit models range from 3.0x10"^ 
Lq for HOPS 108 to l.SxlQ-^ Lq for HOPS 368 & 203 
(see Table HI). The observed Leo (see Table |3] & Fig- 
ure [TU)) computed from the FIR CO lines detected with 
PACS (J = 14 - 13 up to J = 46 - 45) is found to 
be 44% to 81% of the total CO luminosity. Thus, if 
the CO emission from protostars is dominated by hot 
gas (T > 2000-3000 K) at low densities (^(Ha) < lO^-^ 
cm~^), most of the CO luminosity is emitted within the 
PACS wavelength range. The low density, high temper- 
ature solutions obtained from the fits to single compo- 
nent isothermal models do not appear to be driven by 
the high— J CO lines. We modeled rotational diagrams 
for CO lines observed in the Rl spectral band {Jup = 
14—25) separately, and found that, for most sources, a 
single temperature component can explain the observed 
CO emission in the Rl band only for densities n(H2) < 
10^ cm"3 and temperatures T > 2000 K. 

For five sources - HOPS 10, 32, 68, 85 and 343 - where 
fewer than 8 CO lines were detected {Jup < 21), the 
densities are not constrained by the isothermal models. 
These objects are among the lowest Lhoi sources in our 
sample (see Figure [T51) . Although high density, low tem- 
perature LTE solutions cannot be ruled out for these 
sources, they are also consistent with the low density 
(n(H2) < 10^ cm-3), high temperature (T > 2000 K) 
solutions in the context of single component isothermal 



models. 

Two sources in which several high— J CO lines are de- 
tected (HOPS 182 and 370), the highest- J lines (Jup > 
38) are not well fit by the single is othermal comp onent 
models. This difficulty was noted bv lNeufeldl (pOTa ). who 
pointed out that CO rotational diagrams with large pos- 
itive curvatures cannot be reproduced by a single tem- 
perature component model; instead, these can be mod- 
eled using a medium with a power-law distribution of gas 
temperatures. 

5.1.2. Medium with a power-law distribution of gas 
temperatures 

Next, we model the observed CO emission as arising 
from an optically thin medium with a continuous distri- 
bution of gas temperatures to account for the presence of 
an admixture of gas tempe ratures a l ong th e line of sight. 
Following the approach of iNeufeldl ()2012D . the gas tem- 
perature distribution is approximated as a power-law in 
these models: the column density of the medium in the 
temperature interval between T and T -\- dT has the 
form dN oc dT, ov er a temperat ure range of 10 K 
to 5000 K (for details see lNeufeldlf2012D . Such power-law 
temperature distribution models have been successful in 
describing the observed rotationa l diagrams of molecular 
hydrogen [N eufeld & Yuan 2008: lYuan fcNeu feld 2011]; 
IGiannini et al.ll2011i ). We computed the synthetic rota- 
tional diagrams for CO emission from such a medium for 
a large range in density n(H2), the power-law index b, 
and the column density parameter A^(CO) (for details 
seel^ufcld 2012). The best fit solutions were obtained 
by minimizing the reduced-x^. 

The power-law temperature models constrain the den- 
sity of the emitting gas only for three sources, HOPS 108, 
182 & 370, where more than 25 CO lines are detected. 
These models are shown in Figure [T7] and the best fit 
values are listed in Table S) For HOPS 108, whose ro- 
tational diagram has only a modest positive curvature, 
the power-law fit is only marginally better (in terms of 
reduced-x^) than the isothermal fit. The rotational di- 
agrams of HOPS 182 and 370 have larger positive cur- 
vature and for these two sources, power-law tempera- 
ture models provide significantly better fits (in terms of 
reduced-x^) than those obtained from a single temper- 
ature component model. For these three sources, the 
density of the emitting medium is tightly constrained 
between 4 x 10^ cm^'^ and 10^ cm^^, indicating that the 
excitation of CO is not thermal. The total CO luminosity 
derived from the best fit power-law solutions (for tran- 
sitions with Jup = 1 to 80) is slightly higher than that 
obtained from the isothermal model fits for HOPS 108, 
182 & 370 (see Table U]). Comparison with the observed 
Leo (Table [3]) shows that 52-72% of the total CO lu- 
minosity in these objects is emitted within the PACS 
wavelength range. 

For all the other sources with fewer number of detected 
CO lines, gas densities are not well constrained. Accept- 
able solutions could be found for high densities (n(H2) > 
10^ cm~^) with steeper power-law indices {b > 3) and for 
low denisities (n(H2) < 10^ cm~^) with b < 2. The low 
density solutions with & ^ 2 approach the isothermal so- 
lutions described in the previous section. The degenera- 
cies in the multi-component temperature fits are primar- 
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Figure 17. (Left) Rotational diagrams of the observed CO emission with the best fits for the power-law temperature model overlaid 
(solid red line). Downward arrows indicate upper limits to the fluxes of the CO lines which are blended with a nearby line. Open circles 
correspond to 3 cr upper limits for the non-detections. (Right) Reduced-x-^ contours for T and n(H2). The star symbol marks the minimum 
value of the reduced-x'^. The contours corresponding to 68.3% (blue), 95.4% (green) and 99.7% (red) confidence levels are shown. 



ily because fewer CO lines are detected in these sources. 
When more lines are detected the models tightly con- 
strain the physical parameters of the emitting medium, 
as in the cases of HOPS 108, 182 & 370. Nevertheless, 
the ranges of density and power-law index found for these 
three sources also appear to be consistent with the ob- 
served CO emission from most of the other sources in our 
sample. This is demonstrated in Figure [THl Power-law 
temperature models with b in the range of 2.0—3.0 and 
densities in the range of 10^ — 10^ cm~'^ can explain the 
observed rotational diagrams of all the other sources for 
which fewer lines are detected. 

5.1.3. Degeneracies m the model fits 



There are significant degeneracies in the model fits. In 
the five sources where fewer than 8 CO lines detected, 
neither models can rule out high density (n(H2) > 10* 
cm~^) solutions where the emitting gas is in LTE. The 
more flexible power-law temperature models tightly con- 
strain density to 10^ — 10^ cm~^ only for the three sources 
where more than 25 lines are detected. For other sources, 
even though the power-law models with densities in the 
range of 10^ — 10^ cm~'^ are consistent with the observed 
emission, we cannot rule out higher density solutions. 
Moreover, the power-law temperature models that we 
use are at uniform density. In reality, the CO emit- 
ting medium is likely to have a range of temperature 
and density. Models with multiple temperature and den- 
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Figure 18. Synthetic rotational diagrams for a medium with a 
power-law gas temperature distribution with b = 3.0 (blue) and 
6=2.0 (red) for densities n(H2) = 10^ cm^'^ (dashed line) and 
^^(112) = 10^ cm~'^ (solid line) overlaid on the observed rotational 
diagrams. Downward arrows indicate upper limits to the fluxes of 
the CO lines which are blended with a nearby line. Open circles 
correspond to 3cr upper limits for the non-detections. 



sity components can possibly reproduce the observed CO 
emission for LTE conditions. 

5.2. Thermal vs sub-thermal excitation 

The 'universality' of CO rotational temperatures (line 
ratios), and the degeneracies of the simple models, sug- 
gest strongly that the CO excitation is largely either ther- 
mal at low temperatures 300 K for the lines observed 
in the Rl band), or sub-thermal at high temperatures 
(> 1000-2000 K). If the CO excitation is thermal, i.e. 
n(H2) > 10® cm~'^, the rotational temperature should ap- 
proach the physical temperature of the gas and is insen- 
sitive to density. This can be seen from Figure [T21 where 
the average rotational temperatures predicted for the CO 
emission from an optically thin, isothermal medium are 
shown for various gas temperatures and densities (also 
sec Ncufcld 2012). At these high densities, the observed 
distributions of Tri and Tb2b indicate a narrow range in 
the gas temperature for both the components (see Fig- 
ure[T9l). This implies that, if the CO is thermally excited, 
the physical temperature of the emitting gas should be 
very similar for protostars whose luminosities differ by 
two orders of magnitude. Additionally, if the CO excita- 
tion is thermal, the observed positive curvature of the ro- 
tational diagrams indicates the presence of multiple tem- 
perature components. The average rotational tempera- 
tures, Tri, Tb2B and Tbsa (see Table[3|), are significantly 
different in each source, which suggests the presence of 
at least 3 temperature components. Moreover, the ob- 
served distribution of the rotational temperatures Tlri 
and TsRi (see Sect ion H7T|) indicate that even the CO lines 
in the Rl spectral band (J„p=14— 25) require at least 
two temperature components if the excitation is thermal. 
Thus for most sources in our sample 4 or more compo- 
nents would be required to explain the observed emission 
if the emitting gas is in LTE. The narrow range in the 
rotational temperatures observed for the protostars in 
our sample then requires that these multiple tempera- 
ture components of the CO emitting gas remain nearly 
identical for sources with Lboi in the range of ^ 2—217 
Lq. Thus, the primary weakness of LTE solutions is that 
they require heating mechanisms that maintain relatively 
constant temperatures in multiple gas components over 
a range of protostellar luminosity spanning two orders of 
magnitude. 

However, if the CO excitation is sub-thermal, gas at 
high temperatures (T > 1000—2000 K) and moderate 
densities (n(H2) < 10^ cm~^) can reproduce the ob- 
served CO emission from protostars. For low densities, 
n(H2) < 10*'^ cm~'^, the CO rotational diagrams are well 
fit by emission from an isothermal medium with tem- 
peratures in the range of 2000—5000 K. One attractive 
feature of this set of solutions is that in the low den- 
sity limit, the rotational temperatures are insen sitive to 
densi ty and depends only on the temperature (jNeufeldl 
[20T1 also see Figure [T9)) . At densities below lO''^'^ cm 
a large range in the gas temperatures results only in a 
narrow range in rotational temperatures. For example, 
emitting gas at temperatures in the range of 500 K to 
5000 K can reproduce the observed range in Tri as can 
be seen from Figure[T91 In this case, it is somewhat easier 
to explain the observed behavior of rotational tempera- 
tures as a function of Ltd. In protostars with higher Lboi 
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the CO emitting gas is possibly hotter, but the result- 
ing rotational temperatures are not very different from 
those produced by cooler CO gas in low Lboi sources. 
Note, however, that the isothermal solutions are only 
illustrative and multiple temperature components must 
be present; sources which show large positive curvature 
in their rotational diagrams (HOPS 182 and 370), the 
high— J {Jup > 38) lines are not well fit by a single tem- 
perature component. 

Alternatively, we can model the temperature of the 
emitting medium as a continuous power-law distribu- 
tion where the parameters are now the power-law in- 
dex and the density of the g has been done suc- 
cessfuUy for the excitation H7 (iNeufeld &: YuanI 120081: 
I Yuan fc NeufeidI [2011 IGiannini et al.l I2011D . The ob- 
served line ratios of all the sources are consistent with 
power-law models with power-law index, b = 2.0—3.0 and 
with densities in the range of 10^ — 10^ cm^^. These den- 
sities are also below the critical densities of the observed 
CO transitions and the excitation is sub-thermal. Even 
when multiple gas components with a power-law distri- 
bution of temperatures are present, a range of excitation 
conditions (ri(H2) and b) can explain the narrow range 
in the observed rotational temperatures as long as the 
excitation is sub-thermal. Additionally, the CO emission 
observed with PACS is dominated by the high temper- 
ature (T > 1000 K) components, even when lower tem- 
perature components are present. For the best fit model 
for HOPS 108 (n(H2) = 6.0x10^ cm^^ and 6=2.4), the 
gas components with T > 1000 K dominates the emission 
even in the Rl spectral band for CO transitions Jup > 18, 
and components with T > 2000 K contributes more than 
50% of the observed flux for lines {Jup > 27 ) observed 
in B2B and B3A spectral bands. 

In reality, a solution intermediate to those obtained 
with the power-law and isothermal models, with a range 
of temperatures in which the lower temperature (T < 
1000 K) components do not contribute significantly to 
the observed emission is likely to be more robust. Thus, 
if the FIR CO emission observed with PACS is domi- 
nated by hot (T > 1000 K), moderate density (n(H2) < 
10^ cm~^) gas present in the vicinity of protostars, the 
observed rotational temperatures would allow for a rela- 
tively wide range in the temperature and density of the 
emitting gas. The ability to find solutions in the sub- 
thermal limit that do not require multiple temperature 
components, each of which must maintain relatively con- 
stant temperatures over a range of protostellar luminos- 
ity, makes the high temperature, sub-thermal excitation 
the more attractive solution for the observed FIR CO 
emission. 

6. PHYSICAL ORIGIN OF FIR CO EMISSION IN 
PROTOSTARS 

Observations of FIR CO lines {Jup ^ 25) from pro- 
tostars prior to Herschel have suggested that the CO 
emission arises fro m molecular gas heated by non- 
dissociativc shocks f Watson et al. 1985; Cecc arelli et al.l 
[1998: Gianninijjt al. 2001; Nisini ct al. 2002^^ Stud- 
ies with Herschel, however, have invoked two indepen- 
dent mechanisms to explain the CO emission observed 
with PACS: (i) UV- heating of the envelope cavity walls 
producing PDRs and (ii) sm all scale shocks along the 
cavity walls (|van Kempen et al... 2010b: .Fich ct al...2010i : 
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Figure 19. Synthetic rotational temperatures, Tn\ (red) and 
Tb2b (grGcn), computed for CO emission from an optically thin, 
isothermal, uniform density medium (dashed contours). The solid 
contours represent the upper and lower bounds for the observed 
values of Tp^i (red) and T^^B (green). The shaded region (purple 
crosses) indicates gas temperatures and densities for which a single 
component model can simultaneously reproduce both the observed 
Tri and Tb2b- 



IVisser et al.l[20ll) . The primary motivation for invoking 
the PDRs has been the narrow (FWHM < 2 km s"^) 
velocity component seen in the spectrally resolved line 
profiles of the low— J {Jup < 10) CO lines observed t o- 
wards low-mass protostars dvan Kernpen et al.l l2009bl [al: 
Ivan Dishoeck et al.l[2009t lYildiz et al.ll2010Ll2012| ).^ce 
PACS lines are spectrally unresolved, it was assumed 
that this narrow component due to UV-heated cavity 
walls also contributes to the line fiux of the higher— J 
{Jup > 14) CO transitions observed with PACS. Both 
the UV-heated and shock-heated gas components were 
found to be necessary to explain the observed CO emis- 
sion in three parti cular protostars modeled in detail by 
IVisser et &? . (2012). In the following, we discuss how well 
these two mechanisms can explain the observed prop- 
erties of the FIR CO emission from a larger sample of 
protostars. We argue that emission from the UV-heated 
PDRs along the cavity walls is unlikely to be the domi- 
nant component of the observed FIR CO emission. Most 
of the observed emission, instead, is likely to originate 
from low density (7i(H2) < 10^ cm^"^) molecular gas, 
heated in outflow shocks, within 2000 AU from the pro- 
tostars. 

6.1. Origin in PDRs along envelope cavity walls 

In PDR models, energy released from the mass ac- 
cretion onto the protostar from the surrounding disk is 
the source of the U V flux which heats the cavity walls 
(jVisser et al.l I2012D . The accretion luminosity of the 
protostars scales with Lboi llKenvo n fc Hartmannlll995t 
lEvans et al.ll2009l ) and therefore Luv is likel y to increase 
with i ncreasing L^oi. In the PDR models of IVisser et al.l 
(|2012l ). most of the UV-excited emission originates at 
densities between 10^ and 10^ cm^"^ and is therefore 
close to LTE. The total CO luminosity in this case is 
proportional to the total number of the CO molecules 
present (or equivalently mass of the CO gas) for a con- 
stant gas temperature. Thus, in protostars with higher 
Lhoi (and therefore higher Luv)i ^ larger amount of gas 
will be heated, in turn resulting in a higher Lqo- The 
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PDR models can explain the observed correlation be- 
tween Leo and Li,oi , but only if the temperatures of the 
gas components remain the same in all the sources. In- 
deed, the various rotational temperatures computed - 
Tlri, Tspji, Tpji and Tb2b ~ for the protostars in our 
sample are found to be within narrow ranges. If the CO 
excitation is thermal, as is the case in PDRs, the average 
rotational temperature computed for a small range in Ej 
should be representative of the physical temperature of 
the gas component contributing to those CO lines. It is 
difficult, however, to obtain similar gas temperatures in 
PDRs along the envelope cavity walls for protostars with 
Lboi ranging over two orders of magnitude. In PDRs, 
gas-grain collisions and the subsequent thermal emission 
from the dust grains is the dominant cooling mechanism 
for molecular gas densities ri(Il2) > 10® cm^"^, and the 

cooling rate is oc ri(H2)^. Since the photo-electric heating 
rate is oc Ljjv '^■(112), the gas temperature roughly scales 
as T oc Ljjv / n {112) ■ The CO gas temperature can re- 
main nearly constant in protostars with widely ranging 
Lboi values only if the envelope density scales linearly 
with Lboi- Figure [21 however, shows that the envelope 
density at 1 AU, pi, is uncorrelated with L{,oi, so it is 
unlikely that the gas temperatures along the cavity walls 
are similar in all the protostars. 

In Section 14.31 we showed that the luminosity of the 
observed CO emission from an outflow lobe ~ 19" (pro- 
jected distance ^ 8000 AU) away from the protostar is 
a factor of two higher than that observed for the on- 
source emission. If the UV-radiation from the protostar 
were responsible for the heating of the gas, then hotter 
and more luminous CO emission would be expected from 
the on-source position, the opposite of what is observed. 
Moreover, the rotational temperature and luminosity of 
the CO emission from the outflow lobe far away from 
the protostar (projected distance ~ 8000 AU) is indis- 
tinguishable from that observed at the on-source posi- 
tions (projected distance < 2000 AU) of all the other 
sources in our sample, pointing to a common origin in 
both cases. The CO excitation at the outflow lobe is 
likely due to heating of the low density (n(H2) < 10® 
cm~^) gas in the outflow shocks and not from the proto- 
stellar accretion-generated PDRs. 

The UV flux generated by protostellar accretion will 
heat up the dense molecular material close to the pro- 
tostar and this material is likely to contribute to the 
CO emission. Indeed, the narrow velocity component 
(FWHM < 2 km s^^) which originates in the UV- heated 
cavity walls is seen in the line profiles of the CO J = 6—5 
and CO J — 7—6 transitions observed towards low-mass 
proto stars (e.g. Ivan Kempen et al.li2609bl lal: lYildiz et al.l 
|2012| ). For these low— J lines, this narrow component 
dominates the line flux. Similar narrow component is 
also seen in the CO J = 10 — 9 line for the pro tostars 
observed with Herschel/RIFI (jYildiz et al.ll2010D . How- 
ever, in most cases the line flux is dominated by emis- 
sion fro m the broacl (FW HM = 25-30 km s"^) com- 
ponent (|Yildiz et al.l I20100 . Moreover, for the spec- 
trally resolved lines, the broad component originating 
in shock heated gas is seen to increase in strength go- 
ing from J^p = 1 to Jur> = 10 (|Yfldiz et al.ll2Q12l l2Q10l: 
Ivan Kempen et aLll2009bl ). Therefore, for the higher ex- 
citation CO {Jup > 14) lines observed with PACS, emis- 



sion from PDRs is unlikely to be the dominant contrib- 
utor to the total flux. 

6.2. Origin in outflow shocks 
6.2.1. Shocks along envelope cavity walls 

Models for CO emission from shocks along envelope 
cavity walls predict most of the emission that contributes 
to the observed CO line fluxes originates in the high 
density (preshock density > 10® cm~'^) gas (jVisser et al.l 
|2012| ). Our modeling here, however, suggests that the 
CO excitation is sub-thermal and that the observed emis- 
sion likely originates in low density gas (n(H2) < 10® 
cm~^). Moreover, we find that the CO emission proper- 
ties observed for the outflow position > 8000 AU away 
from the protostar and for the on-source positions of the 
protostars in our sample are essentially indistinguishable. 
All these results suggest that the shocked gas producing 
CO emission is likely not located at the denser part of 
the cavity walls close to the protostar. 

6.2.2. CO emission from low density shock-heated gas 

The optimal solutions obtained from both the isother- 
mal and power-law temperature models indicate that 
the observed CO emission from protostars in our sam- 
ple arises in low density (n(H2) < 10® cm~^) molec- 
ular gas which has high temperature (T > 2000 K) 
components. Such high temperatures and low densi- 
ties are typical of gas heated by outfiow shocks and 
have been repor ted by previous studies of FIR CO lines 
from protostars (iCeccareUi et al.lll99l iNisini et alll 19991 : 
iGiannini et al.ll200ll ). The low densities inferred for the 
observed CO emission are more compatible with shock- 
heated gas within the envelope cavity along the molecu- 
lar outfiow or at envelope radius > several 100—1000 AU. 
Since the size of the central spaxel, from which the an- 
alyzed spectra are extracted, corresponds to a projected 
radius of 2000 AU from the protostars, the shock emis- 
sion must come from within this distance. Indeed, com- 
pact molecular outfiows on < 1000 AU scales with high 
outflow momentum rates have been obser ved toward pro- 
tost ars (e.g..iSaiitiago-Garcfa et al. 2009: Lee et al.|[2009l : 
ffakahashi fc Holl2012| ). 

Non-dissociative shocks with shock speed Vg 

> 25 km s~^ can heat the gas to temperatures 

> 2000 K ()Kaufman fc Neufeldl 11991 . In general, 
outflows in protostars produce multiple shocks with 
varying velocities and different shock-front geometries. 
If the densities of the post-shock gas are low, e.g., n(H2) 
< lO"*'^ cm~^, then the isothermal model predicts that 
the hottest (T > 2000 K) gas component dominates the 
CO emission in all the transitions observed with PACS. 
Emitting gas at temperatures of 2000—5000 K can also 
produce the narrow range in the observed rotational 
temperatures. For slightly higher densities, e.g., 10^ 
cm~^ < n(H2) < 10® cm~^, multiple gas components 
whose temperature follows a power-law distribution 
with index b ranging from 2—3 can explain the observed 
FIR CO rotational diagrams. Even in this case, the 
hotter (T > 1000 K) gas components produce most 
of the CO emission (J„p > 18) observed with PACS. 
The post-shock gas in bow-shaped C-shoc ks can have 
a power-law t emperat ure d is tribution (Smith et alj 
[1991L INeufeld fc Yuan, [20081: lYuan fc Ncufcld .20111 
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iGiannini et al.l 1201 Several such bow shocks are 
likely to be present along the supersonic molecular 
outflow from protostars even within a projected radius 
of 2000 AU. Classical bow shocks with a parabolic 
shape are exp ected to produce a power-law index , 
h ~ 3.8 (IS^ith fc Bran d 1993 INeufeld et aO [2001 
INeufeld fc YuanI 120081) . Our results, however, indicate 
that h is in the range of 2—3 for the observed FIR 
CO emission from protostars; a similar range for the 
power-law index have also been found for H2 emission 
from shocked gas (Ncufcld ct al. 2009; Yuan & Neufeldl 
l20Tl[) . As pointed' out by Yuan fc Neufeldl (|20Tll ). the 
lower values of b suggest that either the curvature of 
the shock front is smaller than that of a parabola, or a 
mixture of shocks with different shock-front geometries 
ranging from planar to bow is present along the flow. 

The observed outflow momentum rat es {Mv) in pro- 
tosta r s are found to scale with L^oi (jBontemps et alj 
Il996t iTakahashi fc Holl2012D . This correlation, which 
extends over 3 orders of magnitude in Ljjoi, is primar- 
ily drive n by the increase in th e mass loss rate in the 
outflows (jBontemps et al.iri996D . There is also evidence 
that the mass loss rate (M) from young T Tauri stars and 
proto stars scales with the accretion rate ([Hartigan et alJ 
119951 ; Remming et al. in preparation). Since Lf,oi in- 
creases with accretion luminosity in protostars, the mass 
loss rates from outflows is expected to scale with Lboi, 
and therefore the total mechanical luminosity in shocks 
should also scale with Lboi- Higher cooling rates are 
therefore expected for higher L^oi sources, which, as Fig- 
ure [TU] shows, is tracked by Leo- Thus, shock heating 
of a lower density (n(H2) < 10^ cm~^) gas can at least 
qualitatively explain the observed trend between Lqo 
and Lboi as well as the observed lack of dependence of 
Trot on Lfcoi and T^oi- 

7. CONCLUSIONS 

We have analyzed the emission lines due to the rota- 
tional transitions of CO ( J=14— 13 up to J = 46—45) in 
the FIR (57—196 /im) spectra of 21 protostars in Orion, 
obtained with Herschel/P ACS. The observed CO lines 
originate in the warm and hot gas within a projected dis- 
tance of < 2000 AU from the protostars. We searched for 
correlations between the observed CO emission proper- 
ties and protostellar luminosity, evolutionary status, and 
envelope density. We modeled the CO lines to constrain 
the excitation conditions in the emitting gas and identify 
possible heating mechanisms. Our main conclusions are 
listed below. 

1. The total luminosity of the CO lines observed with 
PACS increases with protostellar luminosity over a 
large range (2—217 Lq). 

2. The CO rotational diagrams of protostars obtained 
with PACS show a positive curvature and the ro- 
tational temperature implied by the line ratios in- 
creases with increasing rotational quantum number 
J. A minimum of 3—4 rotational temperature com- 
ponents are required to fit the observed rotational 
diagram in the PACS wavelength range. The ro- 
tational temperature computed from the lowest— J 
lines {Jup = 14—25) observed, Tri, is found to lie 



within a narrow range for protostars while Ltoi val- 
ues range over two orders of magnitude. A simi- 
lar lack of dependence with L^oi is also observed 
for the rotational temperatures determined for the 
higher— J lines. 

3. The observed CO emission properties {Leo, Trot) 
are uncorrelated with evolutionary indicators and 
envelope properties of the protostars such as Ti,oi 
and envelope density. 

4. If the CO emitting gas is thermally excited, i.e. 
in LTE, multiple gas components at different tem- 
peratures are required to reproduce the observed 
emission. The temperatures of these components 
must somehow remain independent of protostellar 
luminosity over two orders of magnitude. 

5. The observed CO emission can also be modeled 
as arising from sub-thermally excited gas with 
high temperature (T > 2000 K) components. An 
isothermal medium at uniform density can repro- 
duce the observed emission for densities, n(H2) < 
lO''-^ cm-3 and temperatures, T > 2000 K. Ob- 
served emission can also be reproduced by a uni- 
form density medium with a power-law tempera- 
ture distribution (T,„,„ = 10 K & Tmax = 5000 K) 
for densities in the range of 10^— 10^ cm~'^. 

6. Emission from PDRs, produced along the UV- 
heated envelope cavity walls, is unlikely to be the 
dominant component of the CO emission observed 
in the PACS wavelength range. We argue that the 
simplest explanation for both the observed corre- 
lation between ico and Lboi and the invariance 
of CO rotational temperatures between protostars 
is that the CO emission is dominated by sub- 
thermally excited, shock-heated gas at high tem- 
peratures (T > 2000 K) and low densities (n(Il2) 
< 10^ cm~^), located within the outflow cavities 
along the molecular outflow or along the cavity 
walls at radii > several 100-1000 AU. 
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APPENDIX 

A. SPECTRAL ENERGY DISTIRBUTIONS OF PROTOSTARS 

The observed SEDs of the protostars in our sample are shown in Figure | A1[ The J, i/, and Kg photometry, 
when available, is from the Two Micron All Sky Survey (2MASS; iSkrutskie et alJ[2006[ ). Broad band fluxes at 3.6, 
4.5, 5.8, 8.0 and 24.0 um were obtained in a joint survey of the Orion A and B molecular clouds by the Infrared 
Array Camera (IRAC; [Fazio ct al. 2004) and Multiband Imaging Photometer (MIPS; Rickc ct al. 2004) aboard Spitzcr 
(jKrvukova et all 120121: IMegeath et al.l I2O120 . Spitzer Infrared Spectrograph (IRS; Ho uck et all 1200411 spectra also 
shown in Figure [Kl\ were obtained a s part of the GO programs 30859 and 50374 (PI: Megeath, S. T.; Poteet et 
al. in preparation: iPoteet et aLll2011| ). For five sources toward the OMC-2 regio n (HOPS 60, 108, 3 68, 369 & 370), 
photometry at 19.7, 31.4 and 37.1 ^m was obtained with SOFIA/FORCAST (| Adams et al.ll2012| ). Images at 70 
and 160 um of all the protostars in ou r sample were obtai ned with Herschel/PACS as part of the HOPS program 
(|Fischer et all 120101: iStanke et all 120101: iFischer et alll2012[ ). Fluxes at these wavelenths were measured within an 
aperture of 9.6" radius at 70 /xm and a 12.8" radius at 160 /um after subtracting the median sky background estimated 
from an annulus extending from the aperture limit to twice that value (Fischer et al. in preparation). In addition, we 
also obtained submillimeter observations of our sample sources at 350 /Ltm with the Submillimeter APEX Bolometer 
Camera (SABOCA) and at 870 /im with the Large APEX Bolometer Camera (LABOCA) on the Atacama Pathfinder 
Experiment (Stanke et al. in preparation). For the submillimeter data, we used the peak flux density per beam 
observed toward the source to construct the SEDs; the FWHMs of the beams were 7.3" at 350 /xm and 19" at 870 /xm. 
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Figure Al. Observed SEDs of the protostars in our sample. The solid circles represent broad band photometry from 2MASS, Spitzer/TRAC 
& MIPS, SOFIA/FORCAST, Herschel/PACS and APEX/SABOCA & LABOCA. The downward arrow indicates upper limits. The solid 
line represent Spitzer/IRS low resolution spectra. 
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Figure Al. continued 
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Figure Al. continued 



